Basal-plane incommensurate phases in hep structures 
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An Ising model with competing interaction is used to study 
the appearance of incommensurate phases in the basal plane 
of an hexagonal closed-packed structure. The calculated 
mean-field phase diagram reveals various Ig-incommensurate 
and lock-in phases. The results are applied to explain the 
basal-plane incommensurate phase in some compounds of the 
A'A"BX4, family, like K2M0OA, KiWOi, RhiWOA, and to 
describe the sequence of high-temperature phase transitions 
in other compounds of this family. 

PACS: 64.60.Cn, 05.50.-(-g, 64.70.Rh 



I. INTRODUCTION 

The microscopic origin of incommensurate phases in fer- 
roelectrics, magnetic materials, binary alloys and other 
related systems is subject of interest since the early 60th 
(for a review see Ref . Q ) . It is now well established that 
incommensurate modulation in the majority of those sys- 
tems is caused by the frustrating competition between 
different interatomic or interspin forces responsible for 
structural or magnetic ordering. 

The first approach proposed for studying frustration- 
induced incommensurate phases was based on the axial 
next-nearest-neighbor Ising (ANNNI) model with com- 
peting uniaxial nearest-neighbor (NN) and ncxt-nearest- 
neighbor (NNN) interactions in which the structural 
units have been described by binary Ising pseudo-spin 
variables Depending on the coupling parameters 

and temperature, this model was shown to exhibit a rich 
diagram of commensurate and incommensurate modu- 
lated phases when the wave vector jumps between ratio- 
nal and irrational values of the reciprocal lattice period, 
a phenomenon known as a devil's staircase. 

A large variety of incommensurate systems is adequately 
described by the ANNNI model and its analogous [^-^. 
The common feature of these models is that the mod- 
ulation occurs in the uniaxial high-symmetry direction. 
There are, however, examples of incommensurate phases 
where wave vector is perpendicular to the high-symmetry 
axis and can occur in more than one equivalent direction. 
The appearance of this kind of incommensurate modu- 
lation was studied for the simple cubic lattice with com- 
peting NN and NNN interactions along the cubic axes 
and their diagonals 0], for the simple hexagonal lattice 



with NN and NNN interactions [|6|-Q and for distorted 
triangular lattice with only NN interaction 

In this paper we employ the mean-field approximation 
to study the Ising model on the hexagonal closed-packed 
(hep) lattice, where frustration is uniquely related to the 
topology of the lattice and is provided by the in-plane 
NN antiferromagnetic interaction. We show, for the first 
time, that this frustration being stabilized by the out- 
plane NN interaction gives rise to incommensurate phases 
with the modulation vector lying in the basal plane. We 
also study the phase diagram of the system when small 
out-plane NNN interaction is included and show that it 
leads to a rich sequence of phase transitions. The used 
mean-field consideration is complementary to the previ- 
ous cluster-approximation studies of the hep Ising model 
|]l^,0 and of the related hexagonal honeycomb-lattice 
Ising model [|o|jl|. 

Our interest to the hep Ising model is provided by 
the transition sequence in several ionic A'A"BX4 com- 
pounds where the orientational ordering of BX4 tetra- 
hedra drives a series of structural phases including in- 
commensurate states 1^ . The variety of transitions can 
be explained on the basis of the hcp-Ising model where 
the orientational states of BX4 tetrahedra are described 
by two discrete Ising variables as proposed by Kurzynski 
and Halawa |l|,|l|. The review of the recent studies of 
the two-spin Kurzynski and Halawa model is given in 
[ p^ . We demonstrate that the rigorous treatment of the 
model explains the experimentally observed basal-plane 
incommensurate phases in some A'A"BX4 compounds 
and several other features not explored in the previous 
studies. 

The paper is organized as follows. In Sect. 2 the struc- 
tural properties of A' A" BX^ family and the motivation 
of the hep Ising model are considered. In Sect. 3 we 
treat the frustrated hep Ising model within the mean- 
field approximation in order to demonstrate the appear- 
ance of basal-plane incommensurate phases. We derive 
the phase diagram of the system and discuss their rela- 
tion with the low-temperature cluster-approximation di- 
agram of [p^|-p^. In Sect. 4 we discuss the application 
of the obtained results to the experimental properties of 
A' A" BXi compounds at high temperatures. 

Exploring the Ising model, we use the formal terminol- 
ogy of magnetic systems. The properties of A' A" BX^ 
compounds are characterized by the electrostatic interac- 
tion of BXa tetrahedra. Consequently, the para-, ferro- 
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and antiferromagnetic terms of the Ising model corre- 
spond to the para- ferro- and antiferroelectric terms for 
A'A"BX4 compounds. 

II. COMMENSURATE AND 
INCOMMENSURATE PHASES IN A'A"BX4 
COMPOUNDS 

The structures and transition sequences in the 
A' A" BXi family were described in detail in the review 
articles Shortly, the A'A"BXi compounds can 

be presented as a set of BXl~ tetrahedral anions and 
A"+ cations that are regularly placed in the sites of a 
hep structure. The A'^ cations form a simple hexagonal 
lattice. 

The order-disorder transitions in these compounds arc 
related to the ordering of degenerate orientations of BX4 
tetrahedra in the A' A" matrix. A possible degree of 
freedom is given by the vertical up/down orientation of 
the tetrahedron apexes and an another one by the pla- 
nar orientation of the tetrahedra. Orientational ordering 
of the tetrahedra breaks the initial hexagonal symmetry 
Pds/mmc and leads to a sequence of structural transi- 
tions when temperature decreases. In the Kurzyhski and 
Halawa model |^3j the vertical and planar orientations 
are described by two coupled binary spin variables, the 
orientational ordering being provided by the interaction 
between neighboring tetrahedra. The ground states of 
the related two-spins Ising Hamiltonian were shown 
to correspond to the experimentally observed variety of 
phases in the A' A" BX4 compounds. There is a com- 
mon hierarchy in the transition sequence of A' A" BX4 
compounds. The vertical orientational ordering occur at 
higher temperatures (typically of 600 — 'dQQK) than the 
planar one (below 600-ftr). 

Characteristic feature of A'A"BX4 family is the ex- 
istence of the incommensurate modulations associated 
with tetrahedra orientation that often appear as inter- 
mediate phases at structural transformations. They can 
be related either to the planar or to the vertical orien- 
tational degrees of freedom of the tetrahedra and occur 
in the low- and in the high-temperature regions, respec- 
tively. 

The incommensurate phases of the first type have been 
observed in a great number of A'A"BX4 compounds 
[^,^. They have the modulation vector directed along 
the pseudo-hexagonal axis. The appearance of this kind 
of incommensurate phase was proposed to be related ei- 
ther to specif 

ic antisymmetric interaction of the planar orientations 
of tetrahedra in the unit cell or to ANNNI like 

interaction of the tetrahedra in neighboring basal planes 

The incommensurate phases of the second type are the 
subject of the present study. The modulation wave vec- 



tor associated with the vertical tetrahedra orientation 
has the incommensurate component directed in the basal 
plane of the hep structure. This phase occurs in alkali 
molybdates and tungstates K2M0O4, K2WO4, Rb2W04 
in the temperature interval of 590 — 750K |pO|-p^ . There 
is no consistent explanation of this kind of incommensu- 
rate phase although several ideas were proposed in p4[ . 
In this work we explain the appearance of the basal-plane 
incommensurate structure as a result of competitive in- 
teractions between the vertical orientations of neighbor- 
ing tetrahedra localized in the same basal plane and in 
the NN basal planes, as shown in Fig. 1. 

In the high-temperature region the planar orientation of 
the BX4 tetrahedra is disordered and the corresponding 
planar spin variable is equal to zero. Thus, in this re- 
gion, Kurzyhski and Halawa model is effectively reduced 
to the one-spin Ising model where the vertical up/down 
tetrahedra orientations are described 

by the spin variable Si = ±1. The following analysis 
will be restricted to this region. 

The interaction between tetrahedra has an electrostatic 
nature. The uniaxial anisotropy of the hep lattice induces 
dipolar moments of BX4 tetrahedra parallel to the c-axis. 
Therefore, the interaction between the NN tetrahedra is 
dominanted either 

by their induced dipolar moments or by intrinsic oc- 
tupolar moments. In both cases the interaction Jin be- 
tween two NN BX4 tetrahedra localized in the same basal 
plane favors opposite vertical orientations and therefore 
has antiferromagnetic nature. In contrast, the sign of the 
NN out-plane interaction Jout depends on structural de- 
tails of the system, like e.g. ratio c/a, effective charges 
etc.. Thus we consider both the cases of the ferro- and 
antiferromagnetic interaction for Jout- 

The mean-field minimization of the free energy given in 
the next section shows that the incommensurate struc- 
ture does exist in a certain region of interaction parame- 
ters Jin and Jout- It appears, however, that the account 
of only NN interaction leads to a degeneracy between dif- 
ferent phases. To remove this degeneracy, we introduce 
the weak out-plane interaction Jq„j between NNN tetra- 
hedra. We show that this interaction results in structures 
that are found experimentally. 

III. MODEL AND RESULTS 

We use the Ising spin variable S'i = ±1 to explore the 
idea about frustration induced basal-plane incommensu- 
rate phases. The Hamiltonian is written as: 

F = - ^ JijSiSj. (1) 

where 
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Jij — '-^ out 



Jin for the NN in-plane sites 
for the NN out-plane sites 
J'^^f for the NNN out-plane sites, (2) 



as shown on Fig. 1. The negative and positive sign of 
interaction constants correspond to ferro- and to anti- 
ferromagnetic interactions. As discussed in the previous 
Section, we assume that Ji„ > 0, and that Jout and Jq„j 
can take both positive and negative values. It is conve- 
nient to introduce the driving dimensionless parameters: 
K, = Jout/ Jin', A = J'gy^f/Jin, the last one being assumed 
to be smaller than one. 

Note the following properties of the Hamiltonian (|^) 
that will be used later: 

(i) When the ground state for some interaction parame- 
ters K, A is known, the ground state for — k, —A is easily 
obtained by inversion of signs of the spins in alternated 
planes. 

(ii) The three-dimensional hep Ising model can be 
mapped onto a two-dimensional honeycomb hexagonal 
lattice Ising model with NN, NNN, and NNNN interac- 
tions Hi. 

The ground state of the system is found by the mini- 
mization of the free energy: 



F = -fcTlnrrexp(-iJ/fcT). 



(3) 



At T = the problem reduces to minimization of the 
Ising energy JijSiSj over all the possible spin 

configurations. The effective procedure for solve this 
problem was developed in ]lO| . Using the mapping hon- 
eycomb — > hep lattice and adopting the results of 
to our variables, we find that six phases: I, II, III, V, 
XII and XIII (in notations of whose structures are 
shown in Fig. 3, occur at T = 0. Their energies per one 
spin are: 

Ei = J„(3 + 3k + 3A) Eii = J„(3 - 3k - 3A) 
Em = Jm(-1 + K - 3A) Ev = Jm(-1 - K + 3A) 

Exil ^ Jin{-1 + K- X) ExiII ^ J,ni-l- K. + X). (4) 

The coexistence lines of these phases (shown in the Fig. 
2 by the dots) are defined by the equilibrium of their en- 
ergies. An important observation given in [l^ is that the 
infinite number of other degenerate ground states exists 
along these lines. These states become stable at finite 
temperature due to the entropy factor. The finit 

e-temperature cluster-approximation study shows that 
the phase sequence can be quite complex (see and 
Refs. therein). 

The low-temperature approach is useful when the values 
of the spins are assumed to take a fixed value of either 
or —1. Close to Tc, the fluctuations of spins are important 
and the absolute value of the average at ^< Si > can be 
substantially smaller than 1. To consider this regime, we 



minimize the free energy (^ within the mean-field ap- 
proximation, the variables ai being considered as varia- 
tional parameters of the model. The standard mean-field 
treatment [|2^ gives the following expression for the free 
energy: 

F = — Jijajdj + fcry~^ / tanh~ ^ s c?s. (5) 
^ ■ Jo 

In the vicinity of the transition, the absolute values of 
(Ti are much smaller than 1. Expanding in Taylor 
series we obtain the Ising-like expression with additional 
nonlinear terms af: 



kT 1 kT 

^ = ^ E + 2 E ■^^^^^'^1 + T2 E 



(6) 



Note that, unlike to the discrete Ising spins Si = ±1, 
the variables cr^ sweep the continuous spectrum between 
-1-1 and —1. Expression (^) with arbitrary coefficients is 
frequently used as a starting phenomcnological functional 
for considering of incommensurate phases in systems with 
competing interactions |p],|6|-p|Jl^] . We use this expression 
as a basic functional that provides the phase diagram of 
the system. The critical temperature of the transition is 
given by the highest value of T where the functional (^) 
firstly becomes unstable with respect to formation of the 
pattern of non-zero ctj. 

Let us calculate the critical temperature and the phase 
diagram below the transition as a function of the inter- 
action parameters k, A. It is convenient to work with 
the Fourier transformed variables dg = '^^(7iexp{iqri). 
Fourier decomposition of (^ gives the following energy 
per one spin: 

1 hj^ 

/ = 9 E(fc^ + 2 J(9))a,a_, + — ^(a,a_,)2 



kT 



— 2^ (Tqa-q(Jq'(^-q', 



(7) 



where 



J{q) = Jin ■ [£'{ni - Uj) + Ke{ni) cos + Xec{2n,) cosq^], (8) 

and £{ni) = cos qni -|- cos qn2 + cos gri3, 
e'ijii — Uj) — cos q(ni — + cos q{n2 ~ n^) -f cos (7(713 — ni). 

We take ni = (0,5/3,0), ^2 = (0/2,-6/6,0), = 
(—a/2, —b/6, 0) as shown in Fig. 1. 
An important assumption was made that expression 
(^ contains no umklapp terms provided by commen- 
surate modulation of the lock-in phases. Actually, the 
only phases that give this contribution are these with 
modulation vector q = (0,0,0), (0,0,27r/c), (0,27r/6,0), 
(0, 27r/&, 27r/c). Their spin configuration aj corresponds 
to the states I, II, III and V of Fig. 3 with equal in-site 
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amplitudes a = \aj\. Free energy of these states will be 
calculated in a more direct way later. 

The finite-g incommensurate structure becomes stable 
when the coefficient kT + 2J{q) is negative. Softening 
of kT + 2J(q) occurs simultaneously in several symmetry 
equivalent points of the g-space. The resulting state is 
provided by the superposition of the corresponding de- 
generate plane waves that can give either a Iq stripe 
phase, or a multi-g, double-periodic phase. Later we will 
show that the Ig-phase is more preferable. The modula- 
tion amplitude is a complex value satisfying cr^ = ctI^. 
Since the functional does not lock its phase we con- 
sider ffq to be real. 

The transition temperature and the modulation vec- 
tor are defined by kTinc — — min^gi, 2 J(q). We found 
that the qc component of the modulation vector q is 
always commensurate with the reciprocal lattice vector 
c* = 2tt/c and takes the values of either or 2tt/c. The 
incommensurate modulation appears in the basal plane 
along of either the a or 6 symmetry directions. Four 
incommensurate phases {qa, 0,0), (ga, 0, 27r/c), {0,qb,0) 
and (0, Qb, 27r/c) are possible. 

The critical temperatures for the phases {qa, 0,0) and 
{qa,0, 27r/c) are given by 



/ J" 



{k-2X)^ 
1±2A 



(9) 



where the upper sign corresponds to q^ = and the lower 
one to qc = 2t:/c. The modulation vector 



2 , , 



(10) 



changes from qa = (lock- in phases I or II) to qa ^ 27t/ a 
(lock-in phases V or III). 

The critical temperatures for the phases {0,qb,0) and 
(0, qb, 2tt/c) are given by 



kT^/J.,n = -4(±4Aa;^ + Ax^ ± (k - 2\)x'^ 



(11) 



-H4(±k - 3)2: =F (k + A) + 1) 
where x = ±((-3 ± 2A) + [(-3 ± 2\)'^ ~ 8X{k t 3)]^^^)/SX. 



The modulation vector 



Qb 



- arccos x 
b 



(12) 



changes from qb = (lock- in phases I or II) to qb ~2n/b 
(lock-in phases V or III). 

The incommensurate phases can exist only in that re- 
gion of parameters k and A when the arguments of arccos 
in (^0|) and in (12) are between +1 and —1. 

Depending on the interaction parameters n and A, all 
four incommensurate phases can be stable. To find the 
regions of their stability one should compare = 
with the critical temperatures of the lock-in phases I, 
II, III, V, XII and XIII. The later is found directly from 



(^ since the necessary summation ^ Jijaia^ was already 
performed when calculating (^. For the free energy of 
the lock-in phases we obtain 



com — 2 (^-^ 



2i?com)o'^ 



kT 

V2' 



(13) 



where Ecom is the energy (^ of corresponding commen- 
surate phase at T = 0. The critical temperatures of 
transitions are given by 



kTr, 



-2Er. 



(14) 



After calculation of the maximal critical temperature 
from T^jj and Tcom, we obtain the resulting phase dia- 
gram as shown in Fig. 2 by solid lines. The symmetry of 
the diagram with respect to change of sign of the both 
interaction parameters k and A follows from the property 
(i) of the Hamiltonian (0). 

To follow the evolution of this diagram when tempera- 
ture decreases one should solve the infinite system of the 
coupled nonlinear variational equations obtained from 
the mean- field functional (|). The rigorous solution of 
this problem is beyond the scope of our study. However, 
qualitative aspects obtained in a more simple way are 
discussed below. 

Note, first, that the stripe region of incommensurate 
phases in Fig. 2 appears at the same place where the in- 
finite degenerate lines of phase transitions are predicted 
by the cluster calculations at T = (doted lines). This 
indicates, similarly to the ANNNI model ||^,^ , the devil's 
staircase behavior of the transition sequence below Ti^c- 
An infinite number of lock-in phases, where the wave 
vector jumps between different rational multiples of the 
reciprocal lattice periods, appear at low temperatures. 
Their phase boundaries converge to the doted lines at 
T = 0. Phases XII and XIII that are stable at low tem- 
peratures can appear from paramagnetic state only via 
intermediate incommensurate phases. The evolution of 
the incommensurate phases when temperature decreases 
is shown qualitatively in Fig. 4 for 

JLy^t — 0. The analogous phase diagram obtained in 
[|l3| , when the possibility of incommensurate phase was 
not considered, is shown in the same Figure by the dashed 
lines. 

The phase diagram of Fig. 4 is obtained within the 
approximation when the nonlinearity of the functional is 
modeled by the quartic terms of the Taylor expansion 
(^ and when the harmonic plane-wave approximation 
for the modulation is used. The incommensurate phase 
is described by a superposition of n harmonic plane waves 
with temperature independent modulation wave vector, 
corresponding to the symmetry equivalent points in the 
g-space and with equal real amplitudes aq. Either stripe 
phase with n = 1 or double-modulated phases with n = 
2 or 3 are possible; in the 3g-case, qi,q2 and q^ form 
an equilateral 120*^ star. The harmonic approximation 
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is as exact as T is closer to Tine- Variation of q and 
contribution of higher harmonics at T < Tine is expected 
to be small and leads to shift of the phase boundaries of 
the commensurate low temperature phases towards the 
lower temerature like it happens in the ANNNI model 

Under these approximations the energy of the incom- 
mensurate state (n) is written as: 



fine = nk(T - Tinc)cr'l + l^^Ta^ + n{n - l)kT(7'^, (15) 



or, after minimization over aq as: 
Ink 



22n-lT 



(16) 



Note that ( p^ ) is minimal when n = 1, i.e., Iq is the 
most stable incommensurate phase. Comparison of (|l^ ) 
at n = 1 with the energy of corresponding commensurate 
phase obtained from (|l3|), by minimization over a 



J com— 4 coni ) •, 



(17) 



gives the phase sequence below Tine- When temperature 
decreases, the incommensurate phase is stable up to the 
temperature 



Tc — ^inc (3 4" V^) (^inc -^com)i 



(18) 



defined by the condition fine — fcom- Below a first or- 
der transition occurs to one of the commensurate phases 
I, II, III, V, XII, or XIII. The region of < 1 where the 
above approximations are valid is placed above the doted 
line on Fig. 4. Below this line one can get only qualitative 
ideas about the behavior of the transition boundaries. 



IV. APPLICATION TO A'A"BXi COMPOUNDS 

We use now the results of the hep Ising model with com- 
peting interaction to study the high-temperature phase 
transitions in the A' A" BX^ compounds provided by the 
vertical orientation of the apexes of the BX^ tetrahe- 
dra. The calculated mean-field phase diagram (Fig. 2) 
reveals various incommensurate and lock-in phases that 
appear just below the transition and at lower tempera- 
tures. Some of the commensurate phases (I, II, III and V) 
were discovered by Kurzyhski and Halawa ||l^. The oc- 
curring commensurate phases, their correspondence with 
notation of jl^ , their lock-in wave vectors and corre- 
sponding symmetry groups are enumerated in Table I. 

Our calculations reveal the following new features of 
the phase diagram. The Iq basal-plane incommensurate 
phases directed either in the a or 6 crystallographic direc- 
tion appear in the model. The direct first order transition 
between phases I and V and between phases II and III 



is possible when NNN interaction is included. The new 
phases XII and XIII can occur at low temperatures. 

Table II presents the high-temperature transition se- 
quences for several typical A' A" BXj^ compounds and 
their correlation with the ratio c/a. (The data were col- 
lected on the basis of P, p^ , p5| , p6[ . ) Only few compounds 
reveal the high-symmetry parent phase PG^/mmc as- 
sociated with dynamically disordered BX4 tetrahedra. 
One can imagine that in the other compounds the phase 
PQ^/mmc is virtually present above the melting point. 

According to their high-temperature transition se- 
quence and to their ratio c/a the A' A"BXi compounds 
can be classified in a following way (see also 14|): 

(i) Located in the middle part of the Table II the non 
underline compounds with the direct transition (real or 
virtual) from the phase PG^/mmc to Pmcn (phase V). 
Only a few compounds of this numerous group are given; 
for other examples see: |Q,|l^,^ 26 1. Notably, for all the 
compounds of this group the ratio c/a varies from 1.27 
to 1.64. 

(ii) Located in the middle part of the Table II the under- 
line compounds where the phase PSml (phase II) occurs. 
Two of them, CS2SO4 and Rb2S0ii^ demonstrate the re- 
constructive transition Pmcn — P3ml. Compounds of 
this group have the similar ratio c/a as compounds of 

(i). 

(iii) Located in the upper part of Table II the al- 
kali molybdates and tungstates K2M0O4, K2WO4, 
Rb2W04 [p0|-p2[ that have the intermediate incom- 
mensurate phase modulated along crystallographical di- 
rection b. The actual transition sequence there is: 
Pmcn — (0, gb, 27r/c) — Pd^/mmc. To our knowledge, 
there are no compounds with a-directed incommensu- 
rate phase, although this phase already appears in the 
NN approximation. This group has the smallest ratios 
c/a ~ 1.23 — 1.25. The other compounds of this group, 
K2ZnCli, K2C0CI4, K2CoBri, have similar ratio c/a 
but the high-temperature phase PQ^/mmc was not re- 
ported. They can be candidates for the basal-plane in- 
commensurate phase if the melting does not precede the 
Pmcn — Inc — Pd^/mmc transition. 

(iiii) Located in the lower part of Table II compounds 
that demonstrate the phase P63 that is a subgroup of 
P63TOC (phase I). Several A' A" BeF^ compounds and 
KLiMoO^, with virtual P63 — PG^/mmc transition be- 
long to this group. Two other compounds, KLiSOi and 
RbLiCrOi, have the sequence P63 — Pmcn — P&^/mmc 
These are particular systems since the com- 
plete vertical ordering of the tetrahedra occurs only in 
the room-temperature hexagonal phase P63. The or- 
thorhombic phase Pmcn is characterized by a partial 
vertical disorder of tetrahedra. These compounds reveal 
the reconstructive transition I-V. They can probably be 
considered as an intermediate case between classes (i) 
and (iii). Compounds of the group (iiii) have the largest 
ratios c/a ~ 1.64 — 1.70. 
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Following Kurzynski and Halawa we assume that 
sign and magnitude of Jout (and therefore of k) depend 
critically on the ratio c/a. The interactions Jin and J^^^ 
are less sensible to variation of c/a. Compounds of the 
groups (i), (iii), (iiii) reproduce the phase diagram of Fig. 
2 if one supposes that A < —0.5 and that Jout changes its 
sign from negative (ferromagnetic) to positive (antiferro- 
magnetic) when hep lattice goes from its expanded along 
c-axis form with c/a > 1.63 to the contracted form with 
c/a < 1.63. It is interesting to observe that the dipolar- 
dipolar interaction between two NN out-plane BX4 tetra- 
hedra changes 

its sign exactly at c/a = 1.63. The value of the modula- 
tion vector qh in molybdates and tungstates is correlated 
with the ratio c/a in the following way: the smaller the 
ratio c/a the more qb deviates from the lock-in value 27r/& 
of the p 

hase Pmcn p2| . This behavior is consistent with our 
calculations ([l2|). Compounds of the group (ii) seem to 
have a positive A and are placed in the upper right corner 
of the diagram. 

It would be interesting to study the evolution of the 
phase sequences in A' A" BX4 compounds with a contin- 
uous variation of the interaction parameters. One can 
achieve this, e.g., by the application of uniaxial pres- 
sure along the hexagonal axes that slightly changes the 
ratio c/a. In particular, one can expect to obtain the 
commensurate - incommensurate transition, a-directed 
incommensurate phases and the phase XII with symme- 
try Pbca (or their subgroups). The classical compound 
K2Se04 from the group (i), having the smallest ratio 
c/a = 1.27 could be a good candidate to achieve a tri- 
critical Lifshitz point and a basal-plane incommensurate 
phase when submitted to an uniaxial pressure. 

Another interesting result could be obtained by the ap- 
plication of an electrical field E along the hexagonal axis 
that breaks the mirror-basal-plane symmetry Si —Si 
in the Hamiltonian (0). This results in the additional 
invariants E^-cri and E^^af in the functional (|^). 
The first one slightly favors the ferroelectric phase I. 
The second one gives the additional third order term 
'^qi<^q2<^q3 in (Q) whcrc thc vectors qi, q2, qs form the 
equilateral triangle (calculations are analogous to those 
in [|j). This would lead to the splitting of the transi- 
tion (0, (7b, 27r/c) — PQ^/mmc onto two transitions, with 
intervention of thc iq modulated incommensurate phase. 

In our consideration an interaction with elastic degrees 
of freadom that is known to be important in ferroelectrics 
has been neglected. We expect that this coupling breakes 
the discontinuity of the transition from P&^/mmc to 
the low temperature phase as it was observed in sev- 
eral A' A" BXi compounds. The corresponding analysis 
is currently in progress. 
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vais and P. Abelard, Phys. Rev. B39, 3361 (1989). TABLE I. The lock-in wave vectors {qa,qb,qc) and space 

symmetry groups of the commensurate phases appearing in 

the hcp-Ising model (see also Fig. 3). The correspondence 
between the notation of Ref."^ used also in this paper and the 
notations used in Ref.** is given. 



Notations 
of 


Notations 
of'' 


Lock-in vector 
{qa,qb,qc) 


Symmetry 
group 


I 


FP 


(0,0,0) 


PGsmc 


II 


AP 


(0,0,27r/c) 


PSml 


III 


CP 


(0,2V&,0) 


Pmnn 


V 


BP 


(0,27r/6, 2tv/c) 


Pmcn 


XII 




(7r/a,0,0) 


Pbca 


XIII 




(7r/a,0,27r/c) 


Pbna 



" T. Kudo and S. Katsura, Prog. Theor. Phys., 56, 435 
(1976). 

^ M. Kurzyhski and M. Halawa, Phys. Rev., B34, 4846 
(1986). 
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TABLE II. The ratio c/a and the high-temperature se- 
quences of transitions for several A'A"BX4 compounds. 
Question mark signifies that either the information about the 
transitions at higher-temperature is not available or the melt- 
ing occurs. The basal-plane incommensurate phase Inc has 
the modulation vector {0, qb,2n /c). 







Transition sequence 


KiZnCU 


1.23 


... — Pmcn—l 


K2WO4 


1.24 


... — Pmcn — Inc — PGs/mmc 


K2C0CI4 


1.24 


... — Pmcn—1 


K2CoBr4 


1.24 


... — Pmcn~l 


K2M0O4 


1.24 


... — Pmcn — Inc — PGs/mmc 


Rb2W04 


1.25 


... — Pmcn — Inc — PGs/mmc 


K2Se04 


1.27 


... — Pmcn — PGs/mmc 


Rb2Mo04 


1.27 


... — Pmcn—? 


KNaS04 


1.29 


... - P3ml-? 


CS2WO4 


1.29 


... — Pmcn— 7 


Rb2Se04 


1.29 


... — Pmcn — P63 /mmc 


K2SO4 


1.30 


... — Pmcn — PGs/mmc 


CS2M0O4 


1.30 


... ~ Pmcn—? 


K2Cr04 


1.30 


... — Pmcn—? 


K2Mn04 


1.30 


... — Pmcn—? 


NaLiBeFi 


1.34 


... — Pmcn—? 


CS2SO4 


1.37 


... — Pmcn — P3ml—? 


Rh2S04 


1.39 


... — Pmcn — P3ml—? 


CsLiBeFi 


1.62 


... — Pmcn—? 


CsLiS04 


1.62 


... — Pmcn—? 


RbLiS04 


1.64 


... — Pmcn—? 


KLiBeF4 


1.64 


...-P63-? 


KLiW04 


1.65 


... — P63 — Cubic 


KLiMo04 


1.67 


... — P63 — Cubic 


TlLiBeF4 


1.68 


... -P63-? 


KLiS04 


1.69 


... — P63 — Pmcn — PQz/mmc 


RbLiBeF4 


1.69 


... -P63-? 


RbLiCr04 


1.70 


... — P63 — Pmcn — PGs/mmc 



Since the ratio c/a for different compounds is given at dif- 
ferent temperatures the error bars are estimated as ±0.02 as 
the typical variation of thermal expansion. For compounds 
having an orthorhombic Pmcn symmetry this ratio was esti- 
mated as c/{ab/\/?>Y^'^ . 
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